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We study experimentally the intermittent progress of the mechanically-induced martensitic trans¬ 
formation in a Cu-Al-Be single crystal through a full-field measurement technique: the grid method. 
We utilize an in-house especially designed gravity-based device, wherein a system controlled by wa¬ 
ter pumps applies a perfectly monotonic uniaxial load through very small force increments. The 
sample exhibits hysteretic superelastic behavior during the forward and reverse cubic-monoclinic 
transformation, produced by the evolution of the strain field of the phase microstructures. The 
in-plane linear strain components are measured on the sample surface during the loading cycle, and 
we characterize for the first time the strain intermittency in a number of ways, showing the emer¬ 
gence of power-law behavior for the strain avalanching over almost six decades of magnitude. We 
also describe the non-stationarity and the asymmetry observed in the forward vs. the reverse trans¬ 
formation. The present experimental approach, which allows for the monitoring of the reversible 
martensitic transformation both locally and globally in the crystal, proves useful and enhances our 
capabilities in the analysis and possible control of transition-related phenomena in shape-memory 
alloys. 
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I. INTRODUCTION 

Shape-memory alloys (SMAs) are active crystalline 
materials with desirable mechanical properties, exploited 
in a number of applications from engineering to medicine, 
as actuators, sensors and dampeners m- In these 
substances, effects such as super- (or pseudo-)elasticity 
and shape memory originate from a reversible marten¬ 
sitic transformation triggered by either stress or temper¬ 
ature. These phenomena are strongly affected by the 
formation of austenite/martensite microstructures, typ¬ 
ically including twins and habit planes in various ar¬ 
rangements which may act as proper avenues for the 
reversible phase change @hh]. A number of theoreti¬ 
cal approaches have been developed for the analysis of 
such morphologies, going back to the classical ‘crystallo¬ 
graphic theory of martensite’ PSHIE], based on consid¬ 
erations of kinematic compatibility between phases and 
variants. This viewpoint has been much developed in 

In parallel, the experimental investigation of SMA mi- 
crostructure has adopted a variety of means, primar¬ 
ily optical and electron microscopy, originating an ex¬ 
tensive literature, see for instance the references in & 

[3 . Developments in the techniques of full-field measure¬ 
ments have recently also enabled the experimental anal¬ 
ysis of the microstructural strain field over extended re¬ 
gions of a sample’s surface. Moire interferometry and 
digital image correlation are typical examples of such 
techniques, which collect information on the response of 
the tested SMA samples by distinguishing the different 
phases through their strain amplitude [24h28] . Infrared 


recent decades 


thermography has also been employed to measure the in¬ 
homogeneous temperature distributions associated with 
the phase change in SMAs, giving relevant information 
on the features of the martensitic transformation [291135] . 

The spatial distribution of the phases and variants over 
a sample can also be deduced through the grid method 
[36] , wherein this information is derived by analyzing the 
deformation of a periodic grid deposited on the sample, 
obtaining the strain levels which are generally different 
among the phases and their variants. A main advantage 
of the grid method is that it offers a good compromise be¬ 
tween strain resolution and spatial resolution [36] . How¬ 
ever, this full-field measurement technique have not yet 
been employed to investigate phase-change intermittency 
in memory materials. This inhomogeneity is a fundamen¬ 
tal feature of the martensitic transformation in SMAs, 
wherein under a smooth thermal or mechanical driving, 
the strain progresses through a set of discrete avalanche- 
like events, corresponding to transitions between neigh¬ 
boring metastable states. Such effects have been pre¬ 
viously studied through the measurement of calorimet¬ 
ric and/or acoustic signals [37H44] . These techniques 
have revealed many aspects of the bursty character of 
the martensitic transformation, with events exhibiting, in 
many cases, power-law size distributions, possibly after 
training [45]. This links reversible SMA martensites to 
other material systems exhibiting intermittent dynamics 
and scale-free behavior, in magnetism, superconductivity, 
brittle fracture, crystal plasticity j46H50| . in turn framing 
them within the background of a wide variety of complex 
systems in nature exhibiting avalanche-mediated behav¬ 
ior, as in turbulence, earthquakes, computer or social 
networks, and financial markets [5Qh52] . 
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While recent investigations have examined quantita¬ 
tively some spatial aspects in the evolution of marten¬ 
sitic transition phenomena pjgmjmg, there is a lack of 
systematic quantitative information on the strain events 
derived from the analysis of evolving strain maps dur¬ 
ing the phase change in SMAs. In the present study we 
used the grid method to perform the first full-field inves¬ 
tigation of the intermittent progress of the microstruc- 
tural strain field in a Cu-Al-Be single crystal across the 
mechanically-induced martensitic transformation. This 
was done through uniaxial force-controlled loading tests, 
during which we have recorded and analyzed the super¬ 
elastic and hysteretic stress-strain behavior of the sam¬ 
ple, as well as surface strain data in the forward and 
reverse transition. 

Inspired by earlier work on the stress-controlled load¬ 
ing of SMAs in [39] (see also [321151] ). an experimental ap¬ 
paratus based on gravity was used, which applied to the 
sample a low-rate perfectly monotonic force, controlled 
by very slow flow water pumps, at almost isothermal con¬ 
ditions. The loading device was designed to be capable of 
very small force increments, enabling us to highlight the 
intermittency in the alloy, through the quantitative anal¬ 
ysis of the strain field features measured over the sample. 
In contrast to conventional testing machines, as used for 
instance in [33 |34], in the present apparatus the crystal 
may freely adjust its orientation in relation to the vertical 
loading direction. A main reason for this was obtaining 
relatively less complex microstructures, developed in the 
presence of minimal extraneous effects such as external 
friction, lattice plastification, or thermal inhomogeneity, 
to enable us to investigate the strain transformation in¬ 
termittency occurring in the crystal in its most elemen¬ 
tary and basic form. 

This paper is organized as follows. Sect. [Il] presents the 
sample under study, the experimental setup, the load¬ 
ing conditions annd the procedure employed for image 
processing to retrieve strain maps. Sect. m presents an 
analysis of the hysteresis and microstructures deduced 
from the strain evolution recorded during a full loading 
test. Sect. m discusses the observed intermittent char¬ 
acter of the phase transformation. Sect. [V] gives some 
conclusive remarks. Further details on the experimen¬ 
tal apparatus and on data acquisition and treatment are 
given in the Appendix. The Supplemental Material [55] 
comprises four video files showing different aspects of the 
experimental results. 


II. EXPERIMENTAL METHODS 

We tested a Cu Ain.4 Beo .5 (wt. %) single crystal sam¬ 
ple, the same used in the loading experiments described 
in |33j[34], to which we refer for additional information on 
the crystallographic properties and orientation of the lat¬ 
tice in the sample (see the Appendix). This SMA, under 
uniaxial tensile loading, undergoes a reversible marten¬ 
sitic transformation, from a cubic austenite to a long- 
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FIG. 1. Schematic representation of (a) the SMA sample with 
the gauge region in gray, and (b) the testing device. See the 
Appendix for details. 


period stacking monoclinic martensite with twelve vari¬ 
ants, which are all compatible with the austenitic phase 
0119 ]. The two-fold axes in the monoclinic variants orig¬ 
inate from the cubic axes in the austenite [20 . 

Fig. [lja) shows the sample geometry. Aluminum tabs 
were clamped in two pairs of jaws by means of screws 
onto the ends of the sample, to avoid damage or sliding 
once positioned in the testing machine. A bidirectional 
grid (pitch p = 0.2 mm) was transferred using the tech¬ 
nique described in [SB] on the monitored gauge region 
(dimensions 17.78 mm x 38 mm respectively along the 
x- and //-directions) covering nearly all the exposed sur¬ 
face of the SMA sample, with thickness 0.94 mm. See 
the Appendix for more details. 

The loading apparatus is schematically shown in 
Fig-Sb). The sample, which before testing was austenitic 
with some possible residual martensite, was suspended at 
its top end through a ball joint which allows for any ro¬ 
tations. A tank was hung through a cable to the jaws 
at the bottom end of the sample. The tank was filled 
with water during the experiments, enabling us to apply 
a gravity-load in the vertical direction. Loading and un¬ 
loading of the tank were controlled by electronic pumps 
producing a constant and slow flux of water. This allowed 
for a careful stress control of both the forward and re¬ 
verse austenite-to-martensite phase change upon loading 
and unloading, respectively. As transformation stresses 
depend on the temperature, we ensured the stability of 
the ambient temperature T am b = 26.8 d= 0.5 °C during 
the entire duration of the test (more than 45 h). Due 
to the very slow driving, the testing conditions were al¬ 
most isothermal. Before using the pumps, a pre-load was 
applied to the sample to reduce the test duration short¬ 
ening the first part of the elastic climb where martensitic 
transformation events are negligible. This gave the refer- 
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FIG. 2. (Color online) Stress-strain curve and snapshots of 
the strain field e yy at various stages of the loading-unloading 
cycle (austenite is in blue, martensite in red). The stress- 
strain hysteresis regards the average strain s yy vs. the average 
stress cf yy . The green dots mark points on the transformation 
plateaus for reference in the figures below. Also shown are 
maps of all the in-plane strain components and of the angle 
co, corresponding to the average strain value e yy — 0.067 on 
loading (top), and to e yy = 0.057 on unloading (bottom). 
Notice that in the e xx map the austenite is in red and the 
martensite in blue. See the supplemental video #1 [55]. 


ence configuration for the ensuing strain maps. The for¬ 
ward and reverse transformation plateaus were spanned 
through a slowly-varying, monotonic loading controlled 
by adjusting the water flow through the pumps. 

In what follows we specifically report on a test in 
which we have observed the whole forward and reverse 
transformations under these conditions: (a) pre-load of 
34.37 MPa; (b) constant loading rate 1.055 MPa/h up to 



stress (MPa) 


FIG. 3. Hysteresis in the evolution of the martensitic volume 
fraction v vs. the average stress a yy during the stress-driven 
forward and reverse transformation in Fig. [2] The fact that 
v does not reach 1 originates from the possible presence of 
both residual martensite at the lowest loading as well as some 
residual austenite at the highest loading. 


57.29 MPa (duration: about 22 h); (c) constant unload¬ 
ing rate —0.915 MPa/h down to 35.95 MPa (duration: 
about 23h 15 min). Due to the inherent feedback loop in 
the control, very slow monotonic rates such as the ones 
above, obtained with the present device, are not possible 
with conventional testing machines. As a comparison, in 
earlier stress-controlled tests, such as in [39] . the rate was 
about 37 MPa/h in the gauge region of the sample. 

We have gathered about 20,000 images of the sample 
over the 45-hour duration of the test, from which we have 
determined the spatial distribution of the strain levels on 
the sample surface by applying suitable grid-image pro¬ 
cessing. This enabled us to obtain four maps from each 
image during the loading-unloading cycle, giving the in¬ 
plane linear strain fields e xx , e yy , e xy , and the rotation- 
angle field l d about the ^-direction (see Fig. [l] for the 
axes orientation). More details on the data acquisition 
and treatment are given in the Appendix. 


III. HYSTERESIS AND MICROSTRUCTURE 

Fig. H presents the hysteretic stress-strain relation 
recorded during the loading-unloading cycle. Further in¬ 
formation is given in the supplemental video #1 m- 
Fig- [2] also shows the evolving strain maps through the 
forward and reverse transformations, where we see that, 
expectedly, e xx < 0 and e yy > 0. Fig.[3]shows the associ¬ 
ated hysteresis in the evolution of the martensitic volume 
fraction z/, defined as the percentage of the sample sur¬ 
face where e yy exceeds the threshold e yy = 0.05, which 
is about half of the maximum local strain value recorded 
during the test. 

We observe in Fig. [2] the development of the strain 
microstructure as the phase transformation progresses 
from austenite (blue region in the e yy map) towards a 
unique variant of martensite (red region in the e yy map), 
which is selected among the twelve possible variants due 
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FIG. 4. (Color online) Strain clustering in the sample dur¬ 
ing the forward transformation from austenite to martensite. 
Snapshots for selected growing values of the average strain 
£ yy are shown. The color bar indicates the fraction of pixels 
with the given strain values, with grid pitch on the e xx and 
£ yy axes equal respectively to 0.005 and 0.007. See also the 
supplemental video #2 [55] . 


to the suitable alignment of its lattice with the imposed 
load. The presence of a single variant in this test is in¬ 
dicated by the very homogeneous values for the field e yy 
in the red martensitic domain of the sample (with aver¬ 
age approaching e yy = 0.09 near the end of the trans¬ 
formation), and by the fact that the shear component 
e xy has constant sign. The occurrence of a single vari¬ 
ant is also suggested by Fig. [4j showing how the material 
evolves between strain values mostly clustering at two 
levels, corresponding to (strained) austenite and to a sin¬ 
gle (strained) variant of martensite. At the intermediate 
loading configurations in Fig. [4]a trail is noticed between 
the pure-phases’ strain values. This is due to the elas¬ 
tic deformation of the phases’ lattices under while the 
load is present, as well as to the fact that any small scale 
phase domains in the sample are below the spatial res¬ 
olution of the apparatus (see the Appendix), and thus 
give rise to averaged strain values in the corresponding 
pixels. The evolution of the material from the austenitic 
energy well to the martensitic one along the upper trans¬ 
formation plateau can also be viewed in the supplemental 
video #2 [55] . 

The £ yy strain maps in Fig. [2] show that the forward 
transformation develops through needle-like bands which 
progressively enlarge their width and whose length even¬ 
tually crosses the entire sample. Their inclination agrees 
with the theoretical estimates [59] based on the compat¬ 
ibility equations giving the habit-planes for this trans¬ 
formation [2Qj 21,. The snapshots in Fig. [2] also confirm 
the existence of relative lattice rotations between regions 
occupied by different phases, with jumps of the angle uj 



FIG. 5. (Color online) Time evolution of the strain component 
e yy along the mid-section AB of the sample. The detail in the 
upper left corner of the figure shows the evolution of £ yy in 
the interval [Oh, 14h] at a magnified strain scale. 


from +0.03 to —0.03 rad across the austenite-martensite 
interfaces, again in good agreement with the theoretically 
computed values [59,. The £ yy -field in Fig. [2] appears 
more complex along the reverse transformation, possi¬ 
bly involving also a second habit plane [59] (see more 
details below, and the supplemental videos [55]). The 
evolution on a loading cycle is thus hysteretic not only 
in terms of the stress-strain relation, but also in the dis¬ 
tinct microstructural morphologies developed in loading 
vs. unloading. This might be a factor contributing to a 
slight mechanical irreversibility in the phase change. 

Fig.0 shows, as an example, the time evolution of the 
strain component e yy along the central section AB of the 
gauge region. Most of the phase transformation occurs 
within the time window [14 h, 19 h] on loading (the up¬ 
per plateau in the stress-strain curve of Fig.[2|, and [38 h, 
43 h] on unloading, with a total stress increase/decrease 
of about 5 MPa on each plateau. This gives the overall 
slant of the hysteresis loop in the stress-strain plane in 
the present test. The upper left corner of Fig. [5] shows 
the evolution of s yy in the interval [Oh, 14 h] at a mag¬ 
nified strain scale. The strain profiles in Fig. [5] show the 
evolution of the phase microstructures in Fig. [2] where 
we also notice that the end pixels A and B and their 
adjacent regions do not undergo significant strain. 

It is interesting to contrast these results with those 
recorded in earlier tests performed on the same Cu-Al-Be 
sample by means of a conventional uniaxial loading ma- 
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FIG. 6. (Color online) Comparison of the strain fields ob¬ 
tained in the present loading test (a), and in a test performed 
on the same sample with a conventional uniaxial loading de¬ 
vice in [34] (b). In both cases the austenite-to-martensite 
phase transformation is nearly complete: blue and red shades 
indicate the austenitic and martensitic regions in the sample. 
E yy and e yy denote respectively the Hencky and linear strain 
components. The faint horizontal lines in the left snapshot 
are artifacts due to slight pitch variations in the grid, not 
completely eliminated by the data treatment. The small dots 
originate from local grid defects. 




chine, which also exploited the grid method [34]. Fig. [6] 
reports the ^-strain maps in the two cases, near the end 
of the transformation when comparable small fractions 
of austenite remain in the sample (blue regions). The 
martensitic (red) regions show the different microstruc¬ 
tures developed in the two tests. As noted above, the 
present loading device produced a homogeneous £ yy -field 
(Fig. i left), while the conventional machine produced 
a martensitic zone with two alternating strain levels 
(Fig. © right), which we interpret as martensite twin¬ 
ning. The development of twin bands is indeed expected 
in the conventional machine, as the latter does not al¬ 
low for rotations of the sample, forcing the material to 
activate more than one martensite variant to accommo¬ 
date the imposed load. Conversely, a single variant can 
accommodate the loads in the present testing device, as 
the latter allows for an overall rotation to accompany 
the transformation strain in the sample. This helped in 
producing less complex microstructures, for this first full- 
field study of strain intermittency in SMAs. 



FIG. 7. Time evolution of the jumps A e yy recorded at each 
consecutive image for the average strain e yy during loading 
(a), and unloading. The reference green dots mark the trans¬ 
formation plateaus. Reported in red is the corresponding time 
evolution of e yy , see also Fig. [2] The corresponding heavy¬ 
tailed probability distributions for P(A£ yy ) are also shown 
(logarithmic binning). 


IV. INTERMITTENCY 

We study here various aspects of the non-continuous 
progress of the observed forward and reverse marten¬ 
sitic transformation. We see in Fig. [2] that the global 
variable e yy displays a rather continuous evolution at 
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FIG. 8. (Color online) (a) Evolution in time of the jumps 
As yy , evaluated at each consecutive image, along the mid¬ 
section of the sample. See also Fig. [5] (b) Time evolution 
of A e yy at two pixels PI and P2. The time interval is [13 h; 
18 h], where most of the phase transformation occurred on 
loading. The reference green dots delimit the transformation 
strain plateaus. 


the time scale of the entire test (the same holds for v 
in Fig. [3|. However, Fig. [7] shows that the evolution 
in time of the increments A e yy obtained by comparing 
each subsequent image taken during the test, produces 
a very spiky signal, evidencing a pronounced strain in- 
termittency in the phase-transforming sample under the 
slowly-varying monotonic load. Fig. [7] shows the corre¬ 
sponding heavy-tailed probability distributions P(As yy ) 
for the A e yy amplitudes over almost two decades (we 
threshold the values of e yy and Ae yy at 10 —6 , see the 
Appendix). Yet another form of hysteresis in the for¬ 
ward vs. the reverse transformation is highlighted in this 
way, as we see from Fig. [ 7 ] that most of the jumps As yy 
are systematically smaller on unloading than on load¬ 
ing. However, the largest events are more numerous in 
the reverse than in the forward transformation. Analo¬ 
gous behavior is observed for the spiky signal given by 
the phase-fraction jumps Av (not shown). Fig. [ 7 ] also 
evidences that the recorded strain intermittency is non¬ 
stationary along the transformation cycle, whereby there 
is more activity in the evolution of Ae yy near the center 
of both the plateaus of the hysteresis as opposed to their 
ends (more details on this below). 

The intermittent evolution of the global variables such 
as [y and) e yy originates from the intermittency in the 
underlying strain field in the sample. The latter is high¬ 
lighted in Fig. [8j which reports the bursty behavior of the 
strain increments Ae yy at some typical locations on the 
surface (the increments again refer to each consecutive 
image). The main panel shows the space-time evolution 
of As yy along the sample mid-section on loading. The 
lower panels plot the temporal evolution of Ae yy at two 



FIG. 9. Log-log plot of the probability distribution P(As yy ) 
of the strain jumps Ae yy (logarithmic binning) computed for 
all pixels in the sample during the forward (squares) and re¬ 
verse (crosses) transformation. 


pixels PI and P2 of the same section. Strain events occur 
throughout the forward and reverse phase change in the 
crystal, with the larger bursts proportionally likely to in¬ 
clude phase transformation in the lattice as opposed to 
purely elastic local deformation. We notice high-intensity 
(yellow-to-red) localized strain bursts occur throughout 
the sample, indicating the local evolution of the phase mi¬ 
crostructures in the crystal. The largest events in Fig. [8] 
are of the order of 2 x 10 -3 , but the scale is limited to 
less than 5 x 10 -4 to better show the evolution of A e yy 
in the intensity interval where most of the strain bursts 
occur for the two considered pixels PI and P2. 

To quantify statistically the degree of pixel-level bursti- 
ness in the strain-field evolution across the sample, we 
consider the amplitudes of Ae yy recorded at all pixels 
along the entire transformation cycle. We threshold the 
values of Ae yy at 4 x 10 -4 , to stay above the noise level 
of the local strain measurements (see the Appendix). It 
results that the pixel values of Ae yy throughout the sam¬ 
ple span about an order of magnitude, with heavy-tailed 
size distributions which are rather abruptly truncated, as 
shown in Fig. [^separately for the forward and the reverse 
phase change. Each single pixel thus gives to the over¬ 
all strain change a contribution that is bounded by the 
transformation strain between the austenitic and marten¬ 
sitic configurations. The large transformation events 
must thus necessarily exhibit some spatial structure, ob¬ 
tained as the sum of many smaller local events. 

We now characterize the long-range aspects of the 
avalanche-mediated phase transformation. To take into 
account all the components of the evolving strain field, 
we consider the norm \Ae\ = (Ae 2 y + Ae xx + 2As^ y ) 1 / 2 
of the strain jumps rather than the sole Ae yy component. 
We then check for any pixel P on which \As\ is, at any 
time, higher than 4 x 10 -4 . Given any such P, we track 
a strain avalanche as the collection of pixels neighbor¬ 
ing P whereon \Ae\ is also above threshold. In this way, 
about 14,000 events are obtained from the strain evo¬ 
lution data collected along the hysteresis cycle. We re¬ 
call that the time interval between images in the present 
tests is almost ten seconds (see the Appendix), which is 
many orders of magnitude larger than the temporal scales 
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FIG. 10. Avalanching activity during the forward transfor¬ 
mation. (a) Time evolution of the number tla of avalanches 
detected per consecutive image (black) and its running sum 
(red). The green dots mark the reference points on the trans¬ 
formation strain plateaus, (b) Log-log plot of the correspond¬ 
ing probability distribution P(n A ) (logarithmic binning). 


of avalanche durations as recorded for instance through 
acoustic emission [45]. Thus the avalanches evidenced 
here largely pertain only to a single image, and can give 
no indication on the actual time dynamics of the trans¬ 
formation, as each avalanche merges a large and variable 
number of bursty microscale events. 
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We first consider the number tia of avalanches detected 
per image in the sample. On loading, we record a total 
of almost 8,000 events, and about 6,000 events in the re¬ 
verse transformation. In both case avalanching produces 
a non-st at ionary spiky signal, as shown for instance in 
Fig. [To](a) for the direct transformation. The correspond¬ 
ing heavy-tailed statistics in Fig.pT^b) evidence the very 
fluctuating character of these co-operative strain events, 
with tia spanning two orders of magnitude. The supple¬ 
mental video 7^3 [55] highlights the strain avalanching in 
the sample on loading, which helps visualize the bursty 
development of the martensite bands mentioned earlier. 
The supplemental video 7^4 [55] shows the corresponding 
evolution of the strain avalanches, with their more com¬ 
plex microstructures, along the reverse transformation. 

For each avalanche we then consider its size S given 
by the total number of involved pixels, and its magnitude 
M given by the integral of |Ae| over the avalanche. The 
quantity M measured for each strain avalanche gives the 
contribution of such event to the intermittent deforma¬ 
tion detected in the sample at each consecutive image 
as the loading was applied, and, in particular, to the 
sample elongation detailed in Fig. [7] The probability 


FIG. 11. Strain avalanching in the martensitic transforma¬ 
tion. (a)-(b) Log-log plot of the probability distributions 
P(M) and P(S) of the strain avalanche magnitude M and 
the avalanche size S during the transformation cycle (logarith¬ 
mic binning), (c) Log-linear plot of the correlation between 
M and the ratio M/S for the recorded strain avalanches. 
(d)-(e) Values of the exponents a respectively for P(M) and 
P(S), determined by the maximum likelihood method as a 
function of the lower cutoff imposed to the data. 


distributions for both S and M, computed over the full 
cycle are shown in Fig. [TT] While none of the heavy¬ 
tailed statistics for strain intermittency considered so far 
(Figs. [7pl)|) seemingly exhibit power-law character, the 
log-log plots for P(S) and P(M) in Fig. 11 [a)-(b) appear 


in remarkable agreement with the emergence of scale- 
free behavior during the transition process, which is well 
documented by the power-law size statistics for acoustic- 
emission events reported in many cases for the marten¬ 
sitic transformation of SMAs [371441. 

We notice that the power-law character of the distri¬ 
butions of Figs. E3 b ) -(c) is affected by a slightly larger 
frequency of events with sizes of the order of S ~ 10 2 , 





























4r 


loading 



A4 mi n 


FIG. 12. Exponents a of the avalanches magnitude distri¬ 
bution P(M) respectively for the forward (a), and for the 
reverse (b) transformation. The exponent values are deter¬ 
mined by the maximum likelihood method as a function of 
the lower cutoff imposed to the data. 


which are mainly due to spurious avalanches due to grid 
defects. Also, Fig.[lTJc) shows there is a positive correla¬ 
tion between S and M, with a super-linear tendency for 
the avalanches with large S to have larger M . This leads 
to an increase (to about six) of the number of decades 
recorded for the avalanche magnitudes AA in Fig. [Tl|a), 
as compared to those exhibited by the avalanche size S 
in Fig. [ITJb). Furthermore, the defect-related events in 
the small plume of Fig. [ll[c) for S ~ 10 2 , widen the 
corresponding range of A4, and, at the corresponding 
scales in Figs.[ll[a)-(b), locally slightly perturb both the 
distributions PjS) and P(M). This affects the stabil¬ 
ity of these distributions’ exponents, which, as we see in 
Figs. [Tljd)-(e), can both evaluated only as being in the 
range 1.5-1.8. This is compatible with the exponent 1.6 
recently proposed for the size statistics of the optically- 
observed microstructural transformation events recorded 
in [43] at the surface of a Ni-Mn-Ga polycrystal under¬ 
going a temperature-driven martensitic transformation. 

A further hysteretic aspect of the strain evolution 
in our material is manifested by the different expo¬ 
nents of the distributions P(A4) recorded for the strain 
avalanches on loading vs. unloading, see Fig. [T2j We no¬ 
tice that the exponent in the reverse transformation, near 
2, is higher than the one on the forward transformation, 
near 1.5. An exponent asymmetry of this type was also 
recorded in mi for the temperature-driven phase change 
in a NiMnGa SMA. Also acoustic emission data indicate 
a forward-reverse asymmetry in some SMA polycristals 

m- 

We make a final brief analysis of the non-stationarity 
of the transformation process, already remarked from 
Figs. [7] and [lO] We consider six consecutive time sub¬ 
intervals of equal duration over the transformation cycle 
(three for the forward and three for the reverse transfor¬ 
mation). The number of events in each interval are indi¬ 
cated in the insets of Figs, (a)-(c). We see that the non- 
stationarity of the process is marked by a much higher 
strain activity in the central portions of both plateaus, 
as opposed to their end intervals. This is seemingly 
analogous to the non-stationary signal recorded in other 
studies of SMA phase change [43], or in the behavior 
of some hysteretic magnetic systems m- The present 




unloading 



FIG. 13. (Colon online) Non-stationarity of the strain in- 
termittency along the transformation cycle, (a) Probabil¬ 
ity distributions Pi(A 4), i = 1,2,3 (logarithmic binning), 
of the avalanche magnitudes AA during three consecutive 
equal-duration time intervals on the /orward-transformation 
plateau. The inset indicates each interval, with, in parenthe¬ 
sis, the corresponding number of events: 1265, 6024, 374. (b) 
Values of the exponent a for each distribution Pi(A4), deter¬ 
mined by the maximum likelihood method as a function of 
the lower cutoff imposed to the data, (c-d) Data analogs of 
(a-b) for three consecutive equal-duration time intervals on 
the plateau for the reverse transformation, with correspond¬ 
ing number of events: 213, 4798, 614. 















































































































9 


data also show there is a higher avalanching activity 
when the alloy is mostly austenitic than when marten¬ 
site is predominant. The non-stationary character of the 
strain intermittency across the transformation is also ev¬ 
idenced through the probability distributions Pi(M ) of 
the avalanche magnitudes M pertaining to each one of 
the six time intervals, see Figs.[l3ja) and (c). Figs.p^b) 
and (d) show that the exponents of the Pi (Ad) do not all 
have the same value, with the two distributions relative 
to the end intervals of the reverse transformation actually 
not possessing power-law character. As a further aspect 
of the asymmetry in the two phase change directions, 
we have from Figs. [7] and [13] that the recorded above¬ 
threshold avalanches not only are more numerous, but 
are also typically larger on loading vs. unloading. 


V. CONCLUSIONS 

We have performed a full-field study, via the grid 
method, of the progress of the reversible cubic-monoclinic 
martensitic transformation in a Cu-Al-Be shape-memory 
single crystal. The phase change was driven by a stress- 
controlled loading device based on gravity capable of 
very small force increments, which applied to the sample 
a low-rate, monotonic load at largely isothermal condi¬ 
tions. We have thus recorded the stress-strain hysteresis 
in the forward and reverse transformation, and of the cor¬ 
responding strain evolution in the microstructure, i.e. in 
the spatial assemblies of the phases/variants in the sam¬ 
ple. This was done by following the evolution in time of 
the surface strain maps in the sample obtained from the 
deformed grid data. 

In general, due to the breaking of symmetry in the 
phase transition, the martensite exhibits multiple vari¬ 
ants, with different orientations with respect to the par¬ 
ent lattice, and the detailed understanding of the mi¬ 
crostructure is of wide interest in mechanics and materi¬ 
als science, as it is a key element influencing the macro¬ 
scopic properties and behavior of a multiphase crystal 
[H[ 22 j. However, the design of the loading device sought 
to obtain for this study the simplest microstructures in a 
SMA single crystal, involving only the parent austenite 
and one variant of martensite. The very slow driving al¬ 
lowed us to highlight for the first time in a quantitative 
way the strain intermittency that marks the martensitic 
transformation in SMAs. 

We have characterized the non-smooth progress of 
the transformation by first studying the intermittency 
for both the global (body-averaged) and local (pixel- 
level) strain measures. We have thereafter identified 
the organization of possibly long-range events, i.e. strain 
avalanches, and studied their statistical properties (see 
also the supplementary videos [55]). Many earlier ex¬ 
periments based on the acoustic emission accompany¬ 
ing the martensitic transformation in SMAs have shown 
some decades of avalanche-size scaling in the austenite- 
martensite dynamics nuns]. The statistics that we have 


now uncovered for our strain events indicate that power- 
law behavior pertains also to a broad distribution of 
strain-avalanche intensities, covering almost six decades 
well into the macroscopic sample scales. These results 
extend and make more quantitatively precise the results 
recently reported in [43] about the possible scale-free be¬ 
havior of the surface optical activity in SMA martensitic 
transformation. We interpret this strain phenomenology 
as deriving from the evolution of the crystalline sub¬ 
stance within a complex energy landscape whose local 
minimizers represent distinct microstructural morpholo¬ 
gies [saiaimiin}®. The strain avalanches evidenced 
here mark the progress of the material through such min¬ 
imizers under the loading [ 15 ] . The associated display 
of bursty scale-free behavior, as we presently observe, 
may confirm the indications about the critical nature of 
the reversible martensitic transformation in SMAs [37 V 
l40l l42l l43l 115 ] , in analogy to what is observed also in the 
plastic flow of many crystals [49, .64, 65] • We have also de¬ 
scribed a number of other aspects of the observed phase 
change, evidencing especially the non-stationarity of the 
material response along the transformation plateaus, and 
the asymmetric behavior characterized in the forward 
vs. the reverse transformation process. 

Further experimental work along the lines developed in 
the present study can be envisaged in various directions. 
This should advance our understanding, analysis, and 
possible control of phase coexistence and its evolution 
under changing external conditions, and may aid in the 
design of materials with targeted properties. 
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APPENDIX 

The online Supplemental Material [55] accompanying 
this paper comprises four video files illustrating a number 
of different aspects of the phase transformation. Here¬ 
after we give some details on the experimental methods. 

Fig. [l] shows schematically the Cu Ain .4 Beo .5 (wt. %) 
single crystal sample and the testing apparatus. The 
sample dimensions are 17.78 x 38 x 0.94 mm 3 respectively 
along the x-, y-, and ^-directions in Fig. [l] the vertical 
distance between the clamped aluminum tabs being 38 
mm. The monitored gauge region has dimensions 17.78 
mm x 38 mm respectively along the x and y axes. 

The rotation matrix giving the orientation of the 
austenitic cubic single crystal of Cu-Al-Be with respect 
to the reference axes was measured by X-ray diffraction, 
and is given by, in the same reference frame of Fig. [l] (see 
also [34] 5: 
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R = 


0.7502 0.0570 0.6588 

-0.6612 0.0721 0.7467 

-0.0049 -0.9958 0.0918 


We measured the in-plane displacement, the linear 
strain components and the local rotation on the gauge 
region of the sample through the grid method. This 
technique derives these quantities from the images of a 
grid captured by a camera as the sample deforms un¬ 
der the loading. Grids are transferred onto the sample 
surface using the procedure described in [58]. The grid 
pattern (square, with pitch p = 0.2 mm) was printed 
on a thin polymeric substrate using a 12,000-dpi high- 
definition printer and then laid on a thin adhesive layer 
(E504 glue from Epotecny, whose white color optimizes 
the visual contrast with the black grid lines. After cur¬ 
ing for about 40 h at 37 °C, the polymeric substrate is 
carefully removed, but the black ink of the grid lines re¬ 
main bonded on the thin adhesive layer and the sample 
becomes ready for testing. The deformation of the grid 
can reasonably be assumed to be the same as that of 
the surface of the sample as the glue layer is very thin 
(some tenths of mm), and has negligible stiffness com¬ 
pared to that of the sample, making this a non-intrusive 
measurement technique. The maximum strain that can 
be sustained by the grid without debonding or cracking 
is about 18%, which can well accommodate any strain 
reached during the martensitic transformation in our al¬ 
loy. 

In our experiment, the images of the deformed grid 
during the tests were captured by a Sensicam QE cam¬ 
era featuring a 12-bit/1040 x 1376 pixel sensor and a 105- 
mm Sigma lens. The camera was secured on a mounting 
plate whose position was adjusted so that the lines of the 
grid were parallel to the pixels of the sensor along both 
the horizontal and vertical axes. An appropriate magni¬ 
fication allowed to clearly distinguish the black lines of 
the grid. This magnification was adjusted so that one 
grid pitch was encoded with 5 pixels. Three movable 
light guides fed by a KL 2500 LCD cold light source were 
used during image acquisition to obtain a nearly-uniform 
lighting of the grid. About 20,000 images of the deform¬ 
ing sample were obtained over the whole duration of the 
test, one about every 8 sec, with 15-min pauses about 
every 100 min for data recording and refilling the water 
reservoirs which feed the pumps. 

As the grids are the superimposition of two perpen¬ 
dicular straight line patterns, each of them can be con¬ 
sidered as a periodic marking, whose change due to the 
deformation during the test is proportional to the sought 
in-plane displacement [67] : 

_ V A , _ P AA 

u x — 27r Uy — 2 ?r 


where (j) x and <j) y are the local phases along x and 7/, 
respectively. In practice, rf> x and <f> y are deduced from 
the grid images by calculating at any pixel the windowed 
Fourier transform (WFT), see [671 , for which a Gaussian 
window [36] was employed in the present study. The 
possible appearance of local fictitious strains induced by 
slight grid printing defects was strongly reduced by using 
a motion compensation technique described in [36]. 


The above measurement and data-treatment methods 
provide strain maps with a good compromise between 
strain resolution and spatial resolution, a key factor for 
distinguishing details in the strain maps. The spatial 
resolution of the technique, i.e. the shortest distance be¬ 
tween two independent measurements, is considered here 
to the equal to the width of the Gaussian envelope used 
in the WFT, which we consider equal to 6 a according 
to the classical 3-cr rule. As we set a = p, which is 
the smallest value that can be chosen for a |)8|, and as 
5 pixels encode one grid pitch p, we obtain a spatial reso¬ 
lution equal to 30 pixels in this study, that is, 1.2 mm on 
the sample. Random noise propagating from the cam¬ 
era sensor is one of the main sources of disturbance in 
the strain maps. Considering only the latter, and tak¬ 
ing into account that the camera was set to average 128 
frames to provide one grid image, a threshold equal to 
4 x 10 -4 for \Ae\ has been estimated [69l ? ], and has been 
used to separate the avalanches within the noise floor in 
Sect. [[Vi The same value has also been considered for 
the strain components e XX: e yy and s xy . Correspond¬ 
ingly, for the average strain components and their vari¬ 
ations, such as s yy , A £ yy , etc., computed over a whole 
strain map, we have a lower threshold, as we rely on 
4 x 10 5 measurements to calculate a single global quan¬ 
tity for each field. The obtained value was then rounded 
up to be equal 1 x 10 -6 , which was used as threshold for 
the analysis in Sect. [XT] The strain-increment maps show 
in some cases the presence of parasitic fringes, proba¬ 
bly originating from a slight moire effect between grid 
and camera sensor. Such quasi-periodic noise, for which 
there is at present no automatic elimination procedure 
available in the case of strain distributions |71i [72] , may 
shade some actual strain-map features, thus spuriously 
increasing the number of apparently independent strain 
avalanches. This undesirable effect constitutes one of the 
limits of the present analysis. We have however heuris- 
tically performed avalanche identification also by using 
thresholds for \Ae\ lower than the value 4 x 10 -4 per¬ 
taining to the sole sensor noise. We checked directly that 
this only minimally affects the statistical results reported 
in Figs. 10-13, confirming their robustness. 
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